 function res=twobc(ya, yb)
res=[  ya(1)
        yb(1)-1];
 end